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Abstract 

Frustrated magnetism plays a central role in the phenomenology of exotic quantum states. How- 
ever, because the magnetic structures of frustrated systems are aperiodic, there has always been 
the problem that they cannot be determined using traditional crystallographic techniques. Here 
we show that the magnetic component of powder neutron scattering data is actually sufficiently 
information-rich to drive magnetic structure solution for frustrated systems, including spin ices, 
spin liquids, and molecular magnets. Consequently, single-crystal samples are not prerequisite for 
detailed characterisation of exotic magnetic states. Our methodology employs ab initio reverse 
Monte Carlo refinement, making informed use of an additional constraint that minimises variance 
in local spin environments. By refining atomistic spin configurations, we obtain at once (i) a mag- 
netic structure "solution" — i.e. the orientation of classical spin vectors — (ii) the spin correlation 
functions, and (iii) the full three-dimensional magnetic scattering pattern. 
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Frustrated spin systems are characterised by the existence of a macroscopically degen- 
erate manifold of magnetic ground states which suppresses long-range order pQ; topical 
examples include spin ices such as Ho 2 Ti 2 7 [2] and quantum spin liquid candidates such 
as herbertsmithite jHlHj- While the absolute arrangement of spins differs amongst the de- 
generate states accessible to a given system — and hence there is no unique or periodic 
ground state structure — the states share well-defined local correlations that distinguish 
their magnetic structure from the random-spin arrangement of classical paramagnets. The 
importance of understanding these correlations lies primarily in determining the microscopic 
origin of exotic quantum phenomena that emerge from or within frustrated spin systems, 
notably including e.g. the evolution of high-temperature superconductivity from spin liquids 
[5] and the ability of spin ices to support magnetic monopoles [6]. 

Spin orientations within ordinary magnets can usually be determined using a combina- 
tion of neutron scattering experiments and crystallographic analysis of the magnetic Bragg 
diffraction pattern. The suppression of long-range spin periodicity in frustrated magnets 
means that no magnetic Bragg scattering is observed. Instead the magnetic contribution 
to the neutron scattering pattern is a smoothly varying function of three-dimensional (3D) 
reciprocal space, whose symmetry is dictated by the atomic lattice and whose modulation 
depends on the strength and nature of the magnetic correlation functions [1]. There is as yet 
no generic methodology for recovering the correlation functions from the observed scatter- 
ing pattern. Historically, the neutron scattering pattern anticipated from a predetermined 
interaction model is calculated and compared with the experimental single crystal neutron 
scattering data (see e.g. |7]). There are two fundamental problems with this approach. 
First, it relies on anticipating the nature of the interactions responsible for local magnetic 
ordering. Second, it is inherently unfeasible for the very many interesting materials for which 
large single crystal samples simply do not exist. 

We address both issues by developing a model-independent method of magnetic structure 
determination that exploits the information content of the magnetic contribution I(Q) to 
one-dimensional (ID) powder neutron scattering data. Drawing inspiration from the largely 
analogous problem of using nuclear total scattering data to determine the positions of atoms 
in disordered materials such as glasses and liquids [8] , our approach is to fit powder diffraction 
data by refining the orientations of a configuration of semi-classical spins. We proceed to 
show that the extent of information loss during spherical averaging of the single- crystal 
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(3D) magnetic scattering function J(Q) is actually relatively minimal, and that the full 3D 
scattering pattern is recoverable for each frustrated system of varying complexity, symmetry 
and dimensionality that we explore. Moreover the corresponding real-space correlations 
can now be probed directly: the configurations are a "solution" of the frustrated magnetic 
structure. 

For a given system, we generate a supercell of the crystallographic unit cell and assign 
random spin vectors to each magnetic centre. This spin configuration forms the starting 
point for a reverse Monte Carlo (RMC) refinement, which we carry out in the usual manner 
[HI [ID]- In particular, by employing starting configurations with random spin arrangements, 
we make no assumption regarding the form of the magnetic structure (hence Ll ab initio 
structure refinement" [H]). Random spin orientation moves are generated, then accepted 
or rejected according to the RMC algorithm, until the best possible fit-to-data is obtained. 
We then calculate spin correlation functions from the RMC configurations and check for 
evidence that all spins are in equivalent environments (i.e. the nearest-neighbour cross- 
correlation functions are unimodal, see Fig. S4). If such evidence is found we continue the 
RMC refinement with an additional penalty term Xy ar that acts to minimise variance in the 
six local spin correlation functions S a S/3 (a G {x,y},(3 G {x, y, z}), 



where the w a p are empirical weighting terms and the last sum is taken over all magnetic 
atoms i. In essence this step corresponds to searching for the simplest configurations that 
agree with experiment. We have shown elsewhere that the approach of minimising local 
structure variance can be highly effective at improving RMC models of disordered systems 
|12j . Finally, we assess the quality of the models obtained by calculating the full 3D scat- 
tering function /(Q) and also, where appropriate, through inspection of the configurations 
themselves for evidence of local spin correlations. A more complete discussion of our partic- 
ular implementation of RMC, together with the derivation of the spin correlation functions 
and their variances, is given as supporting information. 

As a first case study we demonstrate the recovery of the "ice rules" in spin ice, Ho2Ti207 
|13j . a well-studied material in which magnetic Ho 3+ ions are arranged at the vertices of a 
pyrochlore lattice of corner- sharing tetrahedra. At the simplest level, its magnetic structure 
involves the four moments of each H04 tetrahedron arranged such that two point directly 
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FIG. 1: Spin ice, Hc^l^CV (a) Fits to neutron powder diffraction data: simulated data (top) and 
experimental [H] data (bottom). Input data are shown in black, RMC fit in red, and difference 
(RMC— data) in blue, (b) Single-crystal magnetic diffuse scattering, images as follows: the "ideal" 
spin ice model described in the text (top image, left panel); RMC-derived pattern obtained via 
fitting to simulated powder data (top image, right panel); experimental [15j data (bottom image, 
upper left-hand panel); dipolar model [15] (bottom image, lower left-hand panel); and our RMC- 
derived pattern obtained by fitting to neutron powder diffraction data of Ref. [T3] (bottom image, 
right-hand panel), (c) Real-space spin distribution function £ calculated from the RMC configura- 
tions (see SI), illustrating preferential spin alignment along (111) directions, (d) Spin-ice structure 
(top) and a representative fragment of a RMC configuration (bottom) illustrating recovery of the 
spin ice rules. 

towards and two directly away from its centroid: it is the macroscopic number of ways of 
satisfying this local constraint for the entire structure which leads to the observed absence 
of long-range order. The ID and 3D scattering functions I(Q) and /(Q) derived from such 
a model are given in Fig. 1(a) and (b); these illustrate clearly the apparent contrast in 
information content for which single-crystal data are usually preferred. RMC refinement 
of initially-random spin configurations yields a good fit to data [Fig. 1(a)]. Perhaps more 
surprising is that the 3D scattering function /(Q) calculated from an ensemble of refined 
RMC configurations reproduces the expected pattern remarkably well [Fig. 1(b)]. 

Concerned that this agreement is a fortuitous result of having used simulated data, we 
proceeded to use the experimental powder neutron scattering data of Ref. [H] to drive a 
second RMC refinement. The fit obtained is certainly worse [Fig. 1(a)]; indeed, the ex- 
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perimental data were never intended for quantitative fitting, since they were obtained not 
by using polarised neutron scattering techniques but rather as the difference between to- 
tal neutron scattering data collected above and below the spin ice ordering temperature. 
Nonetheless the 3D magnetic scattering pattern calculated from our refined RMC configu- 
rations shows a strong resemblance to the experimental scattering data of Ref. [15] [upper 
left-hand panel of Fig. 1(b)]. Of particular note is that these new RMC configurations re- 
produce the variation in scattering intensity attributed to dipolar interactions [Fig. 1(b), 
lower left-hand panel], that is unaccounted for in simplistic "ice rules" models (as above) 
but known to occur experimentally [13]. We note also that the limited Q-range of the 
powder-averaged data has not affected appreciably the quality of the 3D scattering pattern 
obtained from the RMC configurations. Taken together with the absence of stringent data 
normalisation, this observation suggests that RMC refinements are surprisingly robust to 
the inevitable limitations of experimental data. 

Access to spin configurations capable of reproducing the full /(Q) function allows ques- 
tions to be asked of the local correlations responsible for its modulation in reciprocal space. 
On the simplest possible level, it is straightforward to calculate a spin orientation distri- 
bution function. For our experimental-data-driven RMC refinement of Ho2Ti207 such a 
calculation clearly reflects the preferential alignment along the (111) directions expected for 
spin ice [Fig. 1(c)]. Furthermore, inspection of any of the refined RMC configurations (one 
of which is included as supporting information) betrays the ice rules themselves: we find 
that approximately 95% of Ho 4 tetrahedra obey the "2-in-2-out" rule [Fig. 1(d)]. To the 
best of our knowledge, this represents the first model- independent, data-driven "solution" 
of the spin ice structure. 

We proceeded to test whether our RMC approach enjoyed similar success across a range of 
frustrated magnets of varying topology, form of magnetic interaction, lattice symmetry and 
dimensionality, focussing our choices on systems of long-standing interest and/or particular 
currency within the condensed matter physics community. We have assembled a visual 
summary of our results — for Heisenberg pyrochlores [16] (including the "hexagonal spin 
cluster" spinel ZnCr 2 04 [13, HE]), gadolinium gallium garnet (GGG) [19], the hyperkagome 
Na 4 Ir 3 8 [201 EI], the hexagonal extended kagome YBaCo 4 7 [7], and XY and Ising 2D 
kagome lattices [2"2l |2"B"] - in Fig. 2, with full details of these refinements and a further 
parallel refinement of GGG driven by experimental data [3] given as supporting information 
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[Fig. S6]. In each case we find that the 3D magnetic scattering function is almost entirely 
recoverable from the ID powder-averaged data. Moreover, despite showing only specific high- 
symmetry cuts of reciprocal space in Fig. 2, one has access to the full volume of reciprocal 
space from the RMC configurations. Historically this has always been difficult, even with 
access to large single crystal samples, by virtue of the time-consuming experimental process 
of assembling 3D reciprocal space maps from 2D data sets. 

Local magnetic order is less well defined in each of these particular frustrated magnets 
than is the case for spin ice; consequently, direct inspection of the configurations is not always 
so informative. We do, however, recover the signature of spin planarity in the 2D systems de- 
spite employing 3D Heisenberg spins in the RMC refinement [Fig. S5] . For all seven systems 
it is possible to extract the spin correlation functions, which reflect the (anti)ferromagnetic 
nature and correlation length of spin interactions [Fig. S3]; we note that in principle these 
could perhaps be used further to drive independent estimation of spin interaction parameters 

Extending our methodology to a third distinct class of frustrated spin system, we consider 
as a final case study the single-molecule magnet Mo72Fe3o [25J. The magnetic Fe 3+ ions 
occupy the corners of an icosidodecahedron, forming a network of corner-sharing triangles 
that is a zero-dimensional analogue of the kagome lattice. Its frustrated magnetic structure 
has been described in terms of a "3-sublattice" model [26j [27], in which three groups of 
ten collinear spins point along axes separated by angles of 120° [Fig. 3(a)]; the absolute 
orientation of the vectors is not uniquely determined. In this way the vertices of each 
triangle in the icosidodecahedron are decorated by one spin from each sublattice. Once again 
our RMC refinements, driven by simulated neutron powder diffraction data, consistently 
produced spin configurations that shared the general features of this magnetic order, with 
only the absolute orientation of the three ordering vectors differing between configurations 
[Fig. 3(b)]. 

The use of RMC methods in the solution of non-magnetic disordered structures has a 
chequered history [25] : in this light the apparent success for magnetic structure solution is 
perhaps surprising. For frustrated magnets, however, the problem of structure solution will 
be heavily constrained by knowledge of the positions of magnetic atoms. We also find that 
the inclusion of the local spin invariance term of Eq. plays a crucial role for systems where 
the degree of local order is strong; for example, repeating our spin ice refinements without 
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FIG. 2: Summary of results for frustrated magnets, (a) Heisenberg pyrochlore antiferromagnet; (b) 
"hexagonal spin cluster" spinel ZnC^CU; (c) gadolinium gallium garnet (GGG); (d) hyperkagome 
Na^^Os; (e) extended kagome YBaCo407; (f) kagome XY system; (g) kagome Ising system. The 
left-hand column shows the RMC fit to simulated powder diffraction data (calculation details given 
as SI). The central column shows 2D cuts of the magnetic scattering function I(Q): those labelled 
"Exact" are calculated from the model spin configurations used to generate the simulated powder 
diffraction data; those labelled "RMC" are calculated from the RMC spin configurations obtained 
by fitting to these powder data. Diagrams of the corresponding frustrated lattices are shown in 
the right-hand column. 
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FIG. 3: Mo72Fe3o single-molecule magnet, (a) The previously-described spin structure [261 [27]; 
(b) a typical RMC magnetic structure solution obtained from simulated powder neutron scattering 
data. Red, green and blue vectors indicate spin orientations of the three sublattices. 

this term we obtain configurations that are convincing neither in their reproduction of the 
/(Q) scattering function [Fig. S2] nor in their adherence to the spin ice rules [Table S3]. 
Our results suggest two particularly significant opportunities for further study. First, direct 
investigation of magnetic correlations in systems not available as single crystals — such as the 
quantum spin liquid k-(BEDT-TTF) 2 Cu2(CN) 3 [29] - becomes a viable prospect. Second, 
one can now develop microscopic pictures of spin structures where there exists a signature 
of strongly correlated magnetism, e.g. "spin stripes" in high-temperature superconductors 
[30] . Consequently, the methods developed here offer a valuable new tool with which to 
explore the magnetic structures of frustrated systems associated with quantum phenomena 
of fundamental importance. 

We gratefully acknowledge valuable discussions with D. A. Drabold, D. A. Keen, M. J. 
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